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Abstract 



We examine diffractive proton-proton scattering p p — > p X and photo- and 
electroproduction of p° mesons p — > p° X, where X denotes a proton or a final 
state, into which the proton can go by diffractive dissociation. Using a functional 
integral approach we derive the scattering amplitudes, which are governed by the 
expectation values of light-like Wegner- Wilson loops, which are then evaluated using 
the model of the stochastic vacuum. For the proton, we assume a quark-diquark 
structure. From the scattering amplitudes we calculate total and differential cross 
sections for high centre of mass energy and small momentum transfer and compare 
with experiments. Furthermore we calculate isovector form factors for the proton 
and the pion within the same model. 
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1 Introduction 



In this article we study inelastic diffractive scattering at high centre of mass (cm.) en- 
ergies y/s > 20 GeV and small momentum transfer squared \t\ < 1 GeV 2 . The low 
momentum transfer implies that one has to apply nonperturbative methods to investi- 
gate these processes. 

Our numerical analysis will be focused on two types of reactions: hadron-hadron scatter- 
ing (especially pp —> pX) and photo- and electroproduction of a p° meson 7 < -*- ) p — > p° X. 
The experimentally observed increase of total cross sections for hadronic reactions with 
the cm. energy |]], starting at about y/s = 10 GeV, has been described phenomenolog- 
ically by Donnachie and Landshoff || in the context of Regge theory. The rise of the 
total hadronic cross sections can be described by a pomeron with an intercept slightly 
larger than one [p], @] • For inelastic diffraction, the pomeron-photon analogy was applied 
in || to relate the cross section of these reactions in a quantitative way to the structure 
functions of deep inelastic electron-proton scattering. A number of other methods to deal 
with inelastic diffractive reactions have been proposed (for reviews see 
A description of soft hadronic high energy reactions, starting from a microscopic level, was 
developed in @ where in the case of an abelian gluon model the pomeron properties were 
related to nonperturbative aspects of the vacuum like the gluon condensate introduced 
by Shifman, Vainshtein and Zakharov |§ . These methods were generalised to QCD in |§ . 
The quantity governing the scattering amplitude was found to be a correlation function 
of two light-like Wegner- Wilson loops |i0|, |TT|1 . These correlation functions are evaluated 
using the model of the stochastic vacuum (MSV) |12j as formulated in Minkowski space 
in [10, 11, 13|. This method has been applied to various reactions, for example exclusive 
vector meson production P^T5|, elastic hadron-hadron scattering |P7|] , and photo- and 
electroproduction of pseudoscalar and tensor mesons 0,[19[ . 



In this paper we will extend the model to the description of inelastic diffractive hadron- 
hadron scattering and vector meson photo- and electroproduction. The hadronic scatter- 
ing amplitudes as derived in PHTT|, |13H2l| anc ^ a brief summary of the basic features of 
the MSV in its Minkowskian formulation will be given in the chapters |2| and |3|. Previous 
results for dcx/di in the case of elastic scattering |L7| are also reviewed briefly. To obtain 



results for single diffractive dissociation in proton-proton scattering and p°-photo- and 
electroproduction, the loop-loop correlation functions will be calculated applying two ap- 
proaches. Moreover we use two different models for the diffractive final state X. In this 
work, the constituent configuration of baryons is assumed to be of the quark-diquark type 



for the reasons given in |22| . Then baryons act as colour dipoles like mesons. The numer- 
ical results of all reactions considered here are presented and compared to experimental 
data. Finally we calculate the isovector form factors for the proton and pion within our 
model in chapter |5|. Our conclusions and a summary are given in chapter |6|. 
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2 The scattering amplitudes 



In this section we present our basic formulae for hadron-hadron scattering and photo- and 
electroproduction. We treat exclusive as well as inclusive processes, where in the inclusive 
case only diffractive dissociation of one of the participating hadrons is considered. 



2.1 Hadron-hadron scattering 

Consider the reaction 

hxiPx) + h 2 (P 2 ) -> hx(P 3 ) + X(P 4 ), 



(1) 



where hi and h 2 are hadrons, X is again h 2 or a diffractive excitation of h 2 and denote 
the four- momenta (Fig. [1]). The hadrons h±, h 2 are modeled as quark-antiquark and 



h 2 (P 2 




h(P 3 ) 



A 4T > X(P A ) 



Figure 1: h\ + h 2 — > hi + X 



quark-diquark wave packets for mesons and baryons, respectively. For the wave functions 
we have chosen a Bauer-Stech-Wirbel ansatz |?3J. The diffractive final state X is modeled 
by a gg-pair (or quark-diquark pair) in a colour singlet state. We then use two approaches. 
In the first method we use free plane waves for the quark and antiquark. Integration over 
all allowed values in phase space and the closure relation then yield all possible diffractive 
final states X, where the case of elastic scattering is also included. The second ansatz, 
applied to confirm the results of the first method and to gain additional insight into the 
structure of the calculated differential cross sections, uses the wave functions of a two- 
dimensional harmonic oscillator where the ground state corresponds to hadron h 2 and the 
excited states to the diffractive excitations of h 2 . Since these eigenfunctions form a basis, 
the calculation of cross sections can be performed as follows: first the cross section for 
one specific excited state with definite quantum numbers n, m is calculated and then the 
sum over all excited states is taken to get the inelastic inclusive diffractive cross section. 



In the framework of the model presented in |TT],|20[ we obtain the scattering amplitude 
for reaction (fj) as 



S 



fi 



5 fi + i(2n) A 5(P 3 + P 4 -Pi- P 2 )T fl , 



fi 



(2is) 



d 2 &Te iq T -b T J difi _ 



(2) 
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Here J,as is the diffractive profile function for which we get two expressions depending on 
which description of the diffractive final state is used in the calculation. For the plane 
wave description we obtain 



-iA 4T -yT 



J™ e (b T ,/) = - J d 2 X T dV J' dz ^31 (XT, Z) V^VMl - Z>) 

W+(^b T + (i - *)x r , XT )w_(-^b T + (i - z')yT, yt) - i) g , (3) 

where is the relative transverse momentum between the quark and the antiquark (or 
diquark) of X (see Appendix For the oscillator description we obtain 



J; 



2d osc. 
diff 



(b, 



d 2 x T dV ['dz [ 1 dz , w 31 {x T ,z)X n > m (y T ,z')v 2 (y T ,z') 
Jo Jo 

W+(ib T + (- - z)x. T , x T )W_(-^b T + (- - z')y T , y T ) - 1^. (4) 



Here X n ' m (yr, z') stands for the two-dimensional harmonic oscillator wave function with 
quantum numbers n, m. Inserting in (W) the ground state wave function X ' leads to the 



elastic scattering amplitude as given in [17 



1-z' 




1 - z 



Figure 2: Orientation and extension of the light-like Wegner-Wilson loops in a projection 
into transverse position space 

In (H) and (|j) 1031 (xt, 2) denotes the profile function for the overlap between initial and 
final state of hadron hi for fixed x^ and z. ^(yr, z') defines the initial state wave function 
of /i2- The light-like Wegner-Wilson loops W± are given by 



W± := ±trV(C±) = itr Pexp(-z^ d^G^y] 



(5) 



4 



where P denotes path ordering and C± is the curve consisting of two light-like worldlines 
for the quark and the antiquark and connecting pieces at ±00. In (|3|),(!) x T and yx define 
the extension and orientation in transverse position space of the two loops representing 
the two hadrons hi and h 2 respectively, z (V) parametrises the fraction of the longitudinal 
momentum of hadron h\ (h 2 ) carried by the quark. The impact parameter is given by 
hx, the light-cone barycentres of the loops are then located at |by + (| — z)xy and 
— |t>T + (| — z')yT, respectively (see |14| and Fig. §). As x-axis for the transverse vectors 



2 1 1 V2 

xy, yt and we choose qr- 

The symbol (■ ■ -)g denotes the functional integration which correlates the two loops. 
In ([3]) the loop- loop correlation function is multiplied with the profile function W31, the 
incoming wave function ip 2 and then integrated over all extensions and orientations of the 
loops in transverse space as well as over the longitudinal momentum fraction z of the quark 
of hadron h±, which does not break up diffractively. In the plane wave description z' is part 
of the specification of the final state and thus appears as arguement of J^'g avc (b T , z') in ([|). 
The phase space integral then includes an integration over z' . When using the second 
method involving the two-dimensional oscillator functions to describe the diffractive final 
state, one has to insert the function X n,m on the r.h.s. of @ and to integrate over z' . 
Thus Jj? s osc -(b T ) depends for given oscillator function X n,m on b^ only. 

2.2 Photo- and electroproduction 

For the reactions 

7 + V -> P° + 

7* + P - P° + X, (6) 

where 7 (7*) denotes a real (virtual) photon, all equations derived in the previous section 
remain unchanged but of course W3i(xt, z) in @ now denotes the profile function for the 
overlap between incoming photon and outgoing meson. We use the wave functions for 
photon and p° meson as in [113]. The reaction can again be depicted as shown in Fig. [I], 
with the incoming hadron h±(Pi) replaced by a (virtual) photon and the outgoing hadron 
^•1(^3) by the p° meson. The upper vertex now symbolises the dissociation of the photon 
into a quark-antiquark-pair and the following formation of the vector meson out of these 
partons. Again the momentum transfer is purely transverse and given by q^. As before, 
X denotes all possible diffractive final states into which the hadron h 2 can go. Taking all 
this into account, the scattering amplitude (H),(@) applies to this type of reaction as well. 
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3 Evaluation of the scattering amplitudes 



3.1 The loop- loop correlation function in the model of the stochas- 
tic vacuum 

Now we perform the functional integrals in (0) and @ making use of the MSV. 
A detailed presentation of the MSV can be found in PD|-P^|, where both the original 
formulation in Euclidian space-time and the analytic continuation to Minkowski space- 
time are discussed. Here we will focus on some fundamental properties which are relevant 
for the evaluation of the expectation value of the loop-loop correlation function. 
The starting point for the model is the correlator of two gluon field strength tensors 
at points x± and x 2 , parallel-transported to a common reference point o along the two 
curves C X1 and C X2 : 

(£^G%(o,x 1 -,C x1 )G^\o,x 2 -,C X2 )) g = \5 ab F^ pa (x u x 2 ,o;C xl ,C X2 ) (7) 

The right hand side depends only on the points x\, x 2 and the two curves C Xl , C X2 , 
the common reference point o can be freely shifted along the curve C\ 2 = C X1 + C X2 , 
where C X2 denotes the curve oriented oppositely to C X2 . Due to colour conservation, the 
correlation function is proportional to 5 ab . In the MSV the assumption is made that 
Fpvpa is independent of the choice of the connecting curve C 12 . Then Poincare and parity 
invariance require Fp Vpa to be of the following form: 



Fp Vp a{z) = j(fiV P fiW - g^9v P ) [kD(z 2 ) + (1 - k)D x ( 



dD^z 2 ) 

~\~\Z a Z y (j p p ZpZyCj^a -\- ZpZpQ va Z a Z^C/vp) y\ Kj- 



dz 2 



1 , 2nrt Ml S 



2 



G2 = ^(g'FF) = {0\^G%(0)G^(0)\0), (9) 

where z — X\ —x 2 . Here G 2 is proportional to the gluon condensate (0|G^(0)(j7 aM!y (0) |0), 
D and D\ are invariant functions normalised to 1 at z — 0, -D(O) = -Di(O) = 1, and n 
is a parameter determining the non-abelian character of the correlator. The properties 
of the functions D and D\ are further specified through the assumption of the MSV 
that for space-like separations those functions rapidly fall to zero on a scale given by the 
correlation length a ~ 0.3 fm. The Fourier decomposition of those functions is given by 



A(: J ) = / (^4 e ^(^)- ( 10 ) 



do 
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A suitable ansatz for D and Di is given in JTT | : 

27(2tt) 4 



le) 



(11) 



D(k 2 ) 
Di(k 2 ) 

with the constant A = 8a/37r. The functions of (|T0|) , (|l~i"|) can be compared to lattice 
calculations fl24"|, ^25] for the Euclidian version of the correlator ([^D and from a fit one can 
extract the following ranges for the parameters G 2 , a, k 



(8a) 2 (k 2 -\- 2 

2 27(2vr) 4 i 

3 (8a) 2 (k 2 -X- 2 + ie) 



nG 2 a 4 



K 

a 



0.39 ... 0.41, 
0.80 ... 0.89, 
0.33 . . . 0.37 fm. 



(12) 



The expectation value of one single parallel-transported gluon field strength tensor van- 
ishes due to colour conservation and the fact that the QCD vacuum has no preferred 
direction in colour space: 



(G(i)) G = 0, 

where we have used the abbrevation 



(13) 



(14) 



Furthermore the MSV assumes the vacuum fluctuations of the field strengths to be deter- 
mined by a Gaussian process. This implies that correlators of more than two gluon field 
strengths factorise and thus the process is completely defined by the second moment of 
its distribution. Due to the assumption of a Gaussian process and colour conservation all 
n-point functions with odd n vanish and we are left with 



(G(l) . . . G{2n)) G = (G{h)G{i 2 )) G . . . {G{i 2n ^)G{i 2n )) G . 

all pairings 



(15) 



Now we will make a cumulant expansion [20] for the loop-loop correlation function 



W+W_) g = (>V + (ib T + (i - z)x T , x r )HL(-ib T + (i - z')y T , y T )) G (16) 

in (^|),(Q) and then evaluate our result in the framework of the MSV. 

To expand the correlation function, we proceed as explained in [17]. First the line integrals 

along the closed loops C± are transformed to surface integrals with the help of the non- 



abelian Stokes theorem where, following the authors of [IT], we choose the mantle of a 
double pyramid as the integration surface. The basis surfaces S± of the two pyramids 
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are enclosed by the two loops C± (see Fig. |3|). The common reference point o is chosen 
to be the apex, where both pyramids touch, and P + and P_ are the mantle surfaces 
of the two pyramids, respectively. Following |17j we rewrite the two traces over 3x3 
matrices occuring in (|l6|) after inserting ([|]) as one trace (Tr2) of a matrix acting in the 
9-dimensional tensor product space. With the definition 

r (n v . r \ — I <%/(°> ® !) for xeP + ( A7 \ 

tM ' ' x) -~\ G«„(o,x-C x )(l®??) for xeP. ' {U) 



we can now write (jT6|) as the expectation value of one ordered exponential in the product 
space, where the integration surface is given by the mantle P = P + U P_ of the double 
pyramid: 

(W + W_) G = iTr 2 (pexp f-| J^ v {x) G t ^(o,x; C x )j 
The cumulant expansion of this expression up to the second term reads 



11 11 

W+(-b T + (- - z)xr,x r )>V_(--bT + ~ ^J^Yr)^ 



Tr 2 exp (^- 9 - J^ v {x) J^{x') (P(G t ^(o,x;C x )G tiPa (o,x';C x ,))) ^ 
Tr 2 exp C 2 (b T ,x T ,y T , z,z'), (19) 
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where C 2 is a 9 x 9 matrix invariant unter £7/(3) colour rotations. As shown in [|I7]] this 
finally leads to 

(W+W-) = |e-** x + (20) 



with 



G_2_7T 2 

24 



x(h T , x r , y T , 2, z') = — — {I(r xq , r yg ) + I(r xq , r yq ) - I(r xq , r y g) - I(r xq -, r yq )} , 



a \ 11 \2 a J / \ vr y- x x \\ vr y -r x 
I(t x , r y ) = k-Y v y v x dv<[ J — y — l - } K 2 ' 



2 



A / V A 



j \ r y- vr x\ \ 2 K ( \r„ - /'IV | 



A 



+ (1 - k)ttA 4 ( ^^j 3 ^ (i^-^ ) . (21) 



Here G 2 , A, defined in fl3|),(|TTD and K 2: s are the modified Bessel functions of second 

and third degree. The vectors with i = x,y and j = q,q are those from the coordinate 
origin to the positions of the quarks and antiquarks (or diquarks) in transverse space as 
shown in Fig. |2|. Separating the real and the imaginary part of the above expression (x 
is a real function) gives 

W+(^b r + (i - 2)x T ,x T )W_(-ib T + (~ - z')y T , y T ) - l\ 

^ cos Qx(b T , x T , y T , 2, + ^ cos QxO>r, x T , y T , z, z')j - 1 

sin Qx(br, x T , yr, 2, z')) + i| sin Qx(b T , x T , y T , 2, j • ( 22 ) 

This is the final result for the correlation function of two light-like Wegner- Wilson loops 
in the matrix cumulant method [ I7| . If we assume \x\ *C 1, ( p2[) reduces to 

(W + W- - 1) G = |-i( X (bT,XT,yT,^,^)) 2 } , (23) 

neglecting terms of order x 3 an d higher. This is the result of the traditional expansion 
method fil l. When computing the numerical results for the cross sections we are interested 
in, we will use both ([22|) and ( p3[ ) and compare with experimental data. In the following 
formulae, we will write {. . . } as an abbrevation for either of the expressions (j2~2"D or (|23|) . 
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3.2 Hadron-hadron scattering 



We now have to specify the hadronic wave functions occuring in 
wave functions we make the following ansatz: 



'2z(l-z) 



<Mx T , z) 



For the hadronic 

(24) 
(25) 

(26) 

The wave functions X n,m consist of the eigenfunctions X n,m of a two-dimensional harmonic 



V?i(x T , z) v?j(x T , z). 
Here 1^ is a normalisation factor given by 

^dz2z(l-z) e" (2 -^ )2/2 < 



oscillator [£6 for the yr-dependence and an additional part for the ^'-dependence as 
in (g§: 



X n > m (y T , z') 
X n > m (y T ) 



! 2z'(l - z' 



I hi 



X"' m (y T ), 



1 



^/((n + m)/2)! ((n-m)/2)! 

d 



'Si 



^ Vt m 



2 \2Sl 



y T dy T 



n-\--m 
2 



'SI 

1 1 1 



Vt m 



ZSi y T dy T 



n — m 
2 



P 2 y 2 T /2 

e imG y 



(27) 



where 9 y is the angle between y^ and q^. Now we put everything together, inserting 
the wave functions (f24])-(|2~7j) and the results (|22|),(|23|) for the correlation function of the 
Wegner- Wilson loops in (J3|) and We can simplify our expressions by exploiting sym- 
metry properties of the wave and correlation funtions. The replacements — > -x^ and 
z — > 1 — z, which exchange the quark with the corresponding diquark (or antiquark in 
the case of mesons) in hadron hi, lead to x ~ ► — X ( see Fig- an d (CT))- On the other 
hand these replacements leave the wave functions invariant and thus the integration over 
x-r and z averages out the sinx-terms of (p2|) when inserted in (§) and (|J). For single 
diffractive dissociation and elastic scattering we can therefore replace (|22|) by 



- 1) 



G 



- COS I -Y 

3 \3 A 



3 C ° S U X 



1 



:28i 



In the expansion method (VV+VV_ — 1) G in (^) is already even under x ~~ > ~X- I n 
our model, therefore, the expression for the correlation function is purely real in (|23|) 
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and only the real part of fl22| ) contributes. The T-matrix element is invariant under the 
exchange of hadron hi by its antihadron. Thus we get only C = +1 (pomeron) exchange 
and no C = — 1 (odderon) exchange. An imaginary part of the correlation function and 
C — — 1 exchange terms both non- vanishing after integration with the overlap functions 
could arise from the inclusion of higher cumulant terms in (|T9|). We note that in double 
diffractive reactions a C = — 1 part of the scattering amplitude will already be present in 
our approximation with terms up to the second cumulant in the matrix cumulant method. 

For the description of the final state X in the oscillator method we can use analogous 
arguements. The simultaneous replacements — > —yr and z' — > 1 — z' and subsequent 
integration over and z' lead to the cancellation of contributions with odd m in (Q) 
because of the existence of a factor e v in X n ' m (y<r). Since for these functions odd 
m only occur for odd n, the sum over all excited states in the calculation of cross sec- 
tions can be reduced to a sum over wave functions with even n and the corresponding m's. 



3.3 Photo- and electroproduction of p° mesons 

As mentioned in chapter |2|, the main difference in our ansatz between the T-matrix 
element of hadron-hadron scattering and the photo- and electroproduction of a vector 
meson is in the wave functions used to describe the participating hadrons and the photon. 
In electroproduction we have to consider three polarisation states for both the incoming 
virtual photon and the outgoing vector meson. Since in our model processes with a change 
in polarisation are strongly suppressed (see |14]|), they will be neglected in the following. 
So we only have to deal with photon and meson wave functions with the same polarisation. 
Furthermore the expressions for the two transverse cases are the same, which leaves us 
with just two combinations: 
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with 



e = ^M^+MW, (31) 

, / m 2 / 1\ 2 \ 

/l/t(*) = A/ 'wTV / ^a^)expf-^-/-^--J J . (32) 

Here Q 2 denotes the virtuality of the photon. Further parameters are the proton charge 
e, the mass of the p° meson M p = 768.5 MeV, its decay constant f p = 150.7 MeV and 
its radius parameters and normalisation constants c^l/t and Al/Tj which are different for 
the longitudinal and the transverse cases. For the values of these parameters and how to 
fix them see [Q. 

In model investigations it has been shown that confinement effects can be very well ap- 
proximated and the photon densities can be extended to real photons by a Q 2 -dependent 
light quark mass ]15] 



mm = / - %) °' 22 GeV f " Q2 < « , Ql = 1.05 GeV* . (33) 
1 for Q 2 > Ql 

Of course for photoproduction, that is for Q 2 = 0, the longitudinal overlap function (|29|) 
vanishes, as it must be since real photons have only transverse polarisation. 



4 Total and differential cross sections 



To calculate cross sections for the reactions considered, we have to fix our free parameters, 
namely those of the MSV: G2, a and k; and those of the wave functions, the extension 
parameter Sh t and the width of the longitudinal momentum distribution Zh v The set 
of MSV parameters used in this work has been established in [17] for the case of the 
matrix cumulant method giving ( |2"2"D and in |TJ[ for the expansion method giving 
(see Table ||). 





matrix method 


expansion method 




(529 MeV) 4 


(501 MeV) 4 


a 


0.32 fm 


0.346 fm 




0.74 


0.74 



Table 1: The parameters of the MSV for the matrix cumulant and the expansion methods 



The values given in Table [I] should be considered as effective values extracted from fits to 
high energy scattering data using two different approximate formulae. Thus the differences 
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between the values in the second and third column of the table can be taken as a theoretical 
error estimate. With fixed parameters the model gives energy independent cross sections. 
It has been shown in JIT) that both the energy dependence of the cross section and of the 
slope parameter b of elastic scattering can be well described by energy dependent hadron 
extension parameters ShX s )- m EH ^ was found that in the framework of the matrix 



cumulant method energy dependent extension parameters can even describe the energy 
evolution of the whole differential elastic cross sections da/dt up to \t\ ~ 1 GeV 2 . For the 
matrix cumulant method we adopt the parametrisation from p 



/ s \ °- 034 

S p (s) = 0.700 fm. (34) 

This was obtained by fitting the total cross section as calculated in the model to the 
soft pomeron part of the Donnachie-Landshoff (DL) parametrisation for cr tot 0. For the 
expansion method we have established a similar connection between S p and s: 

/ s \ °- 028 

S p (s) = 0.624 fm. (35) 

At y/s = 23.5 GeV, for instance, we get S p = 0.868 fm and S p = 0.745 fm from (|34]) and 
(|35|), respectively. The width of the longitudinal momentum distribution of the proton 
has been chosen as z v = 0.4 which gives a best fit to the isovector form factor of the 
proton calculated in the framework of our model (see chapter |]). 



4.1 Proton-proton scattering 
4.1.1 Elastic scattering 

With all parameters fixed, we can now perform the numerical calculations for the scat- 
tering amplitude (0)-(|j). We will start with a short reminder of the results for elastic 
scattering obtained in |L7| as these will be needed further on. As mentioned above, elastic 
scattering can be studied with our formulae when replacing the excited state wave func- 
tion X n ' m in (U) by the ground state wave function X 0,0 , which is identical to (24). For 
s 3> Mp the differential cross section is given by 

1 1 2 

dcJ ei = ~r- a y \Tfi\ dt. (36) 

167T S 2 



In |T4] it has been argued that the Gaussian shaped distribution of the longitudinal mo- 



mentum fraction z (z 1 ) can be replaced by a delta-function centred at z = 1/2 (z' = 1/2), 
since the function \ (0), which determines the shape of the correlation function, depends 
only weakly on z (z'). A numerical investigation of the total cross section calculated from 
the optical theorem shows that the resulting difference for cr tot is smaller than 1%. The 
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0.4 0.5 

\t\ [GeV 2 ] \t\ [GeV 2 ] 

Figure 4: The differential elastic cross section dcr e i/dt [mb/GeV 2 ] at y/s = 23.5 GeV 
calculated using the matrix cumulant method (dashed line) and the expansion method 
(solid line) compared to the experimental data from 



profit one draws from this simplification is a much shorter computation time in the nu- 
merical analysis^ In the following we will make use of this simplification if not explicitly 
stated otherwise. 

In Fig. § we compare the results from the matrix cumulant and expansion methods to 
experiment. The first method, i.e. using (^), gives a reasonable description of the data 
f° r 1^1 ~ 1 GeV 2 over many orders of magnitude but underestimates the data at small 
\t\. The expansion method, i.e. using (|2"3"D, gives a better description of the data for 
\t\ < 0.2 GeV 2 but overshoots the data by orders of magnitude for larger \t\. A fit of 
the form dcr e i/dt = Aexpbt to the differential cross section gives b = 13.8 ± 0.4 GeV -2 
for the matrix cumulant method and b = 10.0 ± 0.2 GeV -2 for the expansion method, 
respectively. From a fit to the experimental data p8| we obtain b = 11.6 ± 0.1 GeV -2 . 
These fits have been performed within the range < \t\ < 0.2 GeV 2 , since the description 
of the data over a larger |t|-range would require an additional term oc t 2 in the exponent 
of the fit. If we calculate the integrated elastic cross section at ^/s = 23.5 GeV we obtain 
cr e i = 5.0 mb in the matrix cumulant method and <r e i = 7.3 mb in the expansion method 
compared to an experimental value of a c \ = 6.81 ± 0.19 mb []]]. The fact that the elas- 
tic cross section calculated by the expansion method is closer to the experimental value 
than the one from the matrix cumulant method is easily understood from Fig. f|b. In the 

1 All occuring integrals in the expressions for the cross sections are calculated using the Monte-Carlo 
integration subroutine Vegas |2^j . As a rule of thumb, one has to use at least 10 datapoints per 
integration variable. The substitution thus decreases the computation time by roughly a factor of 100. 
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Figure 5: The integrated elastic cross section as a function of \/s compared to the experi- 
mental data from [[l j 

integral over da/dt only the region \t\ < 0.2 GeV 2 contributes significantly and there the 
expansion method describes the data better. In the region |*| > 0.2 GeV 2 the result from 
the expansion method is bigger than the experimental result, which has as a consequence 
that the resulting integrated cross section is too big. In the following we will restrict 
our considerations to cr tot , cr el and differential cross sections for |£| < 0.2 GeV 2 . For these 
quantities the expansion method gives a satisfactory description of the experimental data 
and we will use this method throughout the rest of our paper. In Fig. |5] we show cr c i cal- 
culated in this way for 10 GeV < ^/s < 10 TeV. The calculation agrees reasonably well 
with the experimental data of Q . The data are as well from pp- as from pp-experiments. 
Because our approach does not include, in Regge terminology, any non-leading trajecto- 
ries, we cannot distinguish between these two reactions and they are described by the 
same scattering amplitude. 

4.1.2 Single diffractive dissociation 

We now turn to inelastic diffractive scattering p p — > p X . Using the plane wave method 
0))©) we calculate the differential diffractive cross section as 

da diS =(2n) 4 ^\T fl \ 2 d 5 V, (37) 
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where 



d 5 ^ = tA^T"^ xd 2 P 4 Td 2 A 4T cb' (38) 

(27Tj y 4s z'(l — z J 

is the 5-dimensional phase space measure for the three particle final state formed by the 
first proton which remains intact and the quark and the diquark which describe the second, 
diffractivly excited proton. As stated above the description of the diffractive final state X 
by a free quark-diquark pair also includes the case of elastic scattering. To obtain the cross 
section <r s d for single diffractive dissociation, we have to subtract the elastic contribution 
and then multiply by 2 to account for the reaction where the first proton breaks up and 
the second stays intact. We then find for the integrated single diffractive cross section as 
a function of \fs the result shown in Fig. |[ Comparing our results to experimental data, 




10 1 10 2 10 3 10 4 

[GeV] 

Figure 6: The integrated single diffractive dissociation cross section as a function of \fs 



one has to keep in mind that the overall normalisation uncertainty of the experiments is of 
(9(10%). Furthermore the derivation of integrated cross sections from experimental data 
involves extrapolations of the measured data at given values of t and £ = M\ / s to regions 
where no data exist. The extrapolations depend on assumptions on the shape of the t- 
distribution and the shape of the ^-distribution. Different experiments make different 
assumptions and thus the resulting integrated cross sections differ from each other. The 
experimental values on the integrated single diffractive dissociation cross section quoted 
here use £ < 0.05 as an upper bound in the mass distribution |3l|, |33l , except for |32| 
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where the range is extended to £ < 0.2. In our calculation of <r s d we integrate over all 
values of £. Because the mass spectrum obtained in our calculation decreases rapidly with 
increasing £ (see Fig. |8]), our numerical result of the integrated cross section is dominated 
by the low mass region and is not sensible to the integration range being £ < 0.05 or 
£ < 0.2. 

In Table |2] we give the ratio R of the single diffractive dissociation cross section to the 
sum of the single diffractive dissociation and the elastic cross sections from our model and 
from different experiments. For ^/s = 546 GeV and 1800 GeV we have used the values of 
<T e i and <7 s d as quoted by the UA4, CDF and E710 experiments. For the ISR energy range 
20 GeV < y/s < 60 GeV a lot of data exist. Since the cross sections do not vary much 
over this energy range, we have fitted both a e \ and a s d as being proportional to a small 
power of \fs and have then calculated R as a function of y/s using these fits. The quoted 
ISR i?-value in Table [| is then evaluated at an intermediate energy of y/s = 38.5 GeV. 



R = 0sd/(0"el + a sd ) 


y/s GeV 


our calc. 


values calc. from exp. 


23.5 
38.5 
62.3 


0.47 
0.47 
0.46 


0.49 ± 0.07 ISR 


29,3^ 


546 


0.45 


0.41 ± 0.02 UA4 [31] 
0.38 ± 0.01 CDF |2j 


1800 


0.44 


0.33 ± 0.05 E710 
0.32 ± 0.01 CDF 


33] 



Table 2: The ratio R of the single diffractive dissociation to the sum of the single diffractive 
dissociation and elastic cross sections from the model and from experiments 



As can be seen our model predicts that the diffractive dissociation cross section grows 
more slowly with increasing energy than the elastic cross section. This is in qualitative 
agreement with experiment, where an even slower rise of o" S( j compared to a e \ is observed. 
The result for the differential cross section of the single diffractive dissociation is shown 
in Fig. 0. The calculated curve describes the slope of the diffractive reaction quite well 
even for larger values of Therefore the agreement with the experiment is reasonably 
good within the |t|-range considered here. This could however be partly accidental. In 
processes where the proton breaks up, the exchange of hard and semihard gluons will play 
an important role. This exchange is not described by our model which is an approximation 
for the infrared behaviour of QCD. We have seen in elastic scattering that the expansion 
method overestimates the cross section for values of \t\ larger than 0.2 GeV 2 , see Fig. |], 
and this could simulate the expected contribution of hard or semihard gluon exchange 
in the diffractive dissociation reactions. We stress however that the fast decrease of the 
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Figure 7: The differential diffractive cross section d<r s d/dt [mb/GeV 2 ] at yfs 
compared to the experimental data from [E9] 



23.5 GeV 



single diffractive dissociation cross section da s d/dt for values \t\ < 0.2 GeV 2 is a firm 
prediction of our model. 

To check the validity of our description of the diffractive final state by a free quark-diquark 
pair using plane waves, we now apply the second method, which describes the diffractive 
final state X through a sum of wave functions of excited states of a two-dimensional 
harmonic oscillator, as explained above in section |2.1| . In this description, the final state 
phase space is two-dimensional as in the case of elastic scattering and the differential cross 
section is given by 

(n,m)^(0,0) 



with Tfi from (g). The sum runs over all even n for the reasons given in sect. |3.2| , the 
associated quantum number m runs over m = —n, — (n — 2), . . . , n — 2, n. The numerical 
analysis shows that both calculations are in very good agreement to each other and that 
summing up the contributions from values of n < 6 already gives ~ 98% of the result 
using plane waves. 

In the following we will consider the mass spectrum d 2 <7 s d/(d£dt) of the single diffractive 
dissociation reaction at y/s = 23.5 GeV for t = —0.0525 GeV 2 , where £ is the squared 
mass of the diffractive final state divided by s. In our ansatz with plane wave final states 
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£ then is given by 

i . = Ml = A! T + (1 ~ z'K + (4Q) 
s z'(l — z')s 

Here m 9 and m ? are the masses of the quark and the diquark which describe the excited 
proton state. To take into account thresholds the mass for the quark has been chosen to 
be 330 MeV and for the diquark 660 MeV so that the sum roughly gives the proton mass. 
Going back to we recognise that we now can no longer replace the Gaussian shaped 



longitudinal momentum distribution in the wave function (24) for the hadron hi which 
breaks up by a delta function centred around 1/2, as we have done in the calculations 
before, because z' determines the value of £ in (f40|). This was different for the calculation 
of dcr s d/d;£, where we performed an integration over the full range of A iT in phase space 
and were not interested in any particular value of £. As a consequence of the introduction 
of quark masses the integration over z' now does not run from to 1, but the integration 
limits are given by 



2-^2^ T \j4-^r + {~^^) ■ (41) 

This ensures that the mass spectrum starts at M\ = M 2 where M p is the proton mass. 
Our plane wave description of the diffractive final state of course also includes elastic 
scattering. To compare with experimental results on diffractive dissociation we have to 
subtract the elastic contribution. To do so we argue as follows: to obtain the elastic con- 
tribution, we integrate d 2 cx s d/(d£dt) over £ from £ = M^/s up to £i. We then determine 
£i in such a way that the integral gives the value of the elastic differential cross section 
dcr c i/dt. We now interprete the mass spectrum as consisting of the elastic part, which lies 
between £ and £i and the dissociation part, which starts at £i. This procedure allows us 
to separate the elastic and the dissociation contributions. 

The result of the calculation is shown in Fig. | for t = —0.0525 GeV 2 and the cm. energy 
\/s = 23.5 GeV together with the data points from |3(|, where we determined £i to 



be given by £i = 1.63GeV 2 /s fa 0.003. The experimental values are smeared out over 
a certain range of values for £ because of the detector mass resolution function. This 
explains also the data for the unphysical negative ^-values. As a consequence the large 
peak of the diffractive mass spectrum is much more pronounced in our calculation and 
the experimental distribution is flatter around that peak. To compare directly with the 
experiment, we would have to fold our results with the mass resolution function of the 
detector used in the experiments PU|, but unfortunately, this resolution function can no 



longer be reconstructed [34]. We note that our model should give reliable results for small 
£. Indeed, for large values of £ the model seems to underestimate the data considerably. 
But for this £ region we expect, for instance, that our purely nonperturbative treatment 
of the scattering must be supplemented by hard gluon radiation which should lead to high 
invariant masses for the diffractively excited state. Furthermore our calculation treats the 
final state as a quark-diquark pair and therefore here no confinement effects are included. 
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4.2 Proton-pion scattering 

We present calculations for the reaction p7r ± — > pX. Of course, the vacuum parameters 
G2, a, k stay the same but we still have to fix the pion extension parameters S n and z n 
in (plf ). Proceeding as in the case of proton-proton scattering we find for the parameters 
S n = 0.60 fm and z w = 0.5 at ^/i = 19.5 GeV. 

Again we will first take a look at elastic scattering. For a cm. energy of \/s = 19.5 GeV 
we find for the integrated elastic cross sections a e \ = 3.12 mb compared to an experimental 
value of cr c i = 3.30 ± 0.11 mb |3"5fl . Fitting our result for the differential cross section by 
doci/di = Aexpbt we find b = 8.7 ±0.3 GeV~ 2 compared to the experimentally measured 
value of b = 7.9 ± 0.2GeV~ 2 for ir + p and b = 8.4 ± 0.1GeV~ 2 for n~ p scattering, 
respectively [[36 . 



Moving on to the reaction where the pion breaks up diffract ively, we calculate <r s d and the 
R- value, which we define as in the case of proton-proton scattering. We find cr s d = 1.99 mb 
and R = 0.39 where the experimental values are <r s d = 1.90 ± 0.2 mb and R = 0.37 ± 



0.03 [37| which is quite good agreement. 



4.3 Photo- and electroproduction of p° mesons 

In this section we present the results for photo- and electroproduction of p° mesons pro- 
ceeding exactly as described above. Since several investigations on the reaction — > 
p°p in the approach discussed here have been published [|1^|,|16[ we concentrate on the 
single diffractive dissociation reaction ^*>p — > p°X, where X does not include the proton. 
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The energy dependence introduced in chapter 4 cannot explain the strong dependence 
which has been observed at HERA for photons with high Q 2 . We therefore do the calcu- 
lations with the parameters fixed for y/s = 23.5 GeV and concentrate on ratios of cross 
sections , which should be more or less energy independent. For photoproduction we find 

g - d( 7 P - ^ = 0.6 , (42) 
M7P P P) 



which agrees within the errors with the experimental value cr s d/°"ci — 0.5^q ^ 6 from 
which was measured in the range 50 GeV 2 < y/s < 100 GeV 2 . In our model the ratio 
of single diffractive dissociation to elastic p° production depends very weakly on Q 2 . For 
the range 7 GeV 2 < Q 2 < 35 GeV 2 we find 

"^fT™ = 0.54 , (43) 

Cel(7 *P -> P°P) 

which also agrees with the experimental value cr s d/c e i = 0.65 ± 0.24 from [03] (here 
60 GeV 2 < y/s < 180 GeV 2 ). The Q 2 dependence of our model calculations can be fitted 
by (J s d{Q 2 ) oc (Q/Qo)~ 5 '° compared with <r s d(Q 2 ) oc (Q / Qo)~ 5 ' 8±1 ' 9 from experiment. The 
polarisation ratio of the outgoing p° mesons has been calculated in the range 7 GeV 2 < 
Q 2 < 35 GeV 2 : 

a ^ ,[ 0.71 . (44) 



C"L,sd + C"T,sd 

This has to be compared with the experimental result o"L, s d/ (o"L,sd + crx, s d) = 0.79 ± 0.15 



5 Form factors in the model 

In this chapter we will study form factors within our model. We do not intend to perform 
a precision calculation of form factors but we will apply the calculation to extract values 
for the width of the longitudinal momentum distributions of the proton and the pion, z p 
and z n , respectively, by fitting our results to experimental data. 

5.1 The electromagnetic form factors of the proton 

The coupling of the electromagnetic current to the proton can be described by 

(p(P',s')\j»(0)\p(P,s))=u s ,(P') 



Fi P (Q 2 )Y 



2M n 



■F 2p (Q 2 



(45) 
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where the momentum transfer is q = P' — P, Q 2 



-q 2 , M p is the proton mass and 



F\pi F%p are the Dirac and Pauli form factor of the proton, respectively. Now we choose a 
coordinate system such that q is purely transverse: 



-P+< + -P^ 
2 + + 2 

Ip,< + -P_n M 



2 y 
1 , 





q T 




/ i \ 

o 

o 



^ + M p 2 )/P + 



In the high energy limit, P + — > oo, we get for the matrix element 



(p(P',s')\f(0)\p(P,s)) = P + n^xl, 



Fi P (Q 2 



a qT • cr 

2M„ 



P 2p (Q 2 



(46) 

0) 

X„ + 0(1), (47) 



sec 



where Xsi Xs' are the Pauli two-component spinors. F\ p multiplies the spin-non-flip part, 
F 2p the spin-flip part of the matrix element. Calculating the spin average of this expression 
leads to 



~ $>(P', s)\f(x)\p(P, s)) = P + n» + F lp {Q 2 ) + 0(1). 



(46 



We describe the calculation of the Dirac form factor of the proton in the framework of 
our model in Appendix |B[ In the following we consider the matrix element of the third 
component of the isospin current j'3. Its matrix element between proton states is as 
in (|4"5|) , (|4"7|) , with Fi P replaced by P„, related to the form factors of proton and neutron 
by 



F iv = ~ {F IP (Q 2 ) - F in (Q 2 )) 



1,2). 



(49) 



With the wave functions f^4|) we obtain 



F lv (Q 2 
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— I dz2z(l-z)e- (z -^ 2 / 2z te-^ n2 



(50) 



p Jo 



where I p is the normalisation factor (|26|). For this calculation we need only the expectation 
value of one Wegner- Wilson loop. The expectation value over one single loop is 1 in both 
the matrix cumulant method and the expansion method. Thus, in our model the form 
factor is just the Fourier transform of the squared wave function. 
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We will now use fl5DD to determine z p and S p . It turns out that in the range < Q < 
0.5 GeV the form factor depends sensitively on S p but only weakly on z p . From a fit to 



experiment in this region we obtain S p = 0.77 fm. 



v 

With S p fixed to this value we show 
in Fig. |9| our result (|50"D for F\ v for different values of z p . The experimental values have 
been calculated from the experimental data for Ge p and Gm p from |42[[43| and a fit of 
the experimental data on Geu and Gmu [HI according to (f49|) and the relation between 
the Dirac (Fi p ^ n ) and the electric {GE P , n ) and magnetic (GM P , n ) form factor of the proton 
and neutron, respectively: 



F 1N {Q 2 



G E n{Q 2 ) + tG mn {Q< 
1 + r 



Q 2 



4M 2 N 



(N = p,n). 



(51) 



z p , which fixes the width 



The best fit is found for z p = 0.4. As can be seen from Fig. 
of the longitudinal momentum distribution of the constituents, plays no important role 
f° r Q ^ 0.5 GeV. For larger values of Q however, our fit is substantially improved when 
using a Gaussian shaped ^-dependence instead of a delta-function centred around z = 1/2, 
which is equivalent to z p — > 0. 




Figure 9: The isovector form factor of the proton for S p 
z p compared to the experimental data from |j4~2"HI~4[ 



2 

Q [GeV] 

0.77 fm and different values of 



It has to be noted, that the proton extension parameter S p obtained from (|50| ) is not, 
and need not be, the same as the one used in the hadronic scattering processes in the 
previous chapters. Whereas the hadronic extension parameter has been allowed to be 
energy dependent (see (j34l),(|35"l)) to account for the rise of er tot with ^/s, the extension 
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parameter connected with the form factor has a fixed value for all energies as the form 
factor itself is energy independent and is related to the electromagnetic radius of the 
proton as follows. Using the definitions 



2\p 



-6 



dG Ep {Q 2 



dQ 2 



Q 2 =0 



(52) 



relations (|49[) , (pT| ) and the experimental value 

dG En (Q 2 ) 



dQ 2 



0.019 fm 2 



(53) 



from thermal-neutron-electron scattering [45], we get from our model 



HL = 0.81 fm. 



(54) 



This coincides with the value one obtains for the proton electromagnetic radius when 
describing the electric form factor of the proton by the dipole parametrisation j|2 , which 
also results in r^ m = 0.81 fm. From scattering experiments one finds r^ m = 0.88±0.03 fm or 
r em = 0.92±0.03 fm, depending on which fit is used for the experimental data on Ge p {Q 2 ) 
for small Q 2 f42j. The Lamb shift measurements |^6| give r^ m = 0.890 ± 0.014 fm. Thus 
our result (54), as well as the one calculated from the dipole parametrisation, is smaller 
than the experimental value for r^ m . Our calculation as well as the dipole fit describe the 
data jKJ for Ge p rather well for Q > 0.4 GeV. But for smaller Q the data indicate a 
rapid change in the slope dG e p {Q 2 ) / dQ 2 which is described neither by our model nor by 
the dipole parametrisation. Such an "anomalous" behaviour of Ge p and Ge u for small 
Q 2 has been related to QCD vacuum effects in [47 . 



5.2 The electromagnetic form factor of the pion 

For the charged pions 7r ± the matrix element of the electromagnetic and the third com- 
ponent of the isospin current are equal. Choosing again the coordinate system as in (|46|) 
with M p replaced by m, we get 

(n + (P')\f(0)\n+(P)) = (P + < + P^F^Q 2 ). (55) 

Here the matrix element can be expressed by only one form factor F n . The calculation of 
this matrix element in our model leads to 

F n (Q 2 ) = f dz2z{l-z)e~ {z - 1 i )2/2z *e-^ slQ2 . (56) 

■Mr JO 
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Figure 10: The pion form factor for S n = 0.68 fm and different values of z v compared to 
the experimental data from [H 



We compare (|56|) to experimental data for F n from [48 in Fig. |10|. As for the proton the 
transverse extension parameter S n can be fitted in the range < Q < 0.5 GeV with the 
result S n = 0.68 fm. Using the analogue of relation ( [521) f° r the pion, this value gives 
an electromagnetic radius r^ m = 0.64 fm, which is consistent with the experimental value 
r em = 0.663 ± 0.006 fm J48|. For values Q > 0.5 GeV our fit becomes sensitive to the 
width of the longitudinal momentum distribution of the constituents. For the pion, the 
best fit for the width of this distribution is given by z n = 0.5. The broader distribution 
compared to the proton is related to the smaller mass of the pion, which is in agreement 
with the parametrisation of the hadron wave functions in 



6 Conclusions 

In this work we have calculated total and differential cross sections for elastic and in- 
elastic diffractive scattering at high cm energies and small momentum transfer. In our 
model we start from a microscopic description of the scattering of quark- ant iquark and 
quark-diquark wave packets and use functional integral methods to get expressions for 
the scattering amplitudes. The correlation functions of light-like Wegner- Wilson loops 
governing these amplitudes are evaluated in the framework of the model of the stochastic 
]. The hadron-hadron scattering amplitudes are obtained by multiplying 



vacuum 



10 



the parton scattering amplitudes with suitable hadronic wave functions p3| , or the photon 
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wave function WM in the case of photo- and electroproduction. Both a matrix cumulant 



expansion for the correlation function of two Wegner- Wilson loops as developed in [17 



and an expansion method ||TT|,|T4| are used. 

The free parameters of our model are those of the model of the stochastic vacuum: G2, 
a and n, and the ones in the wave functions: Sh f and Zh determining the width of the 
transverse and longitudinal momentum distributions of the constituents of the hadrons, 
respectively. These parameters have been determined in previous work [11, 17] on elastic 
scattering. The extension parameters are allowed to depend on the cm. energy 
according to fl34| ) and (^) respectively. In this sense different hadrons are characterised 
through their radii, which come out close to the corresponding electromagnetic radii of the 
hadrons for energies ~ 20 GeV. The values for are obtained from a calculation of 
form factors in our model. Our result for the isovector Dirac form factor of the proton and 
the electromagnetic form factor of the pion, as well as the electromagnetic radii extracted 
from them, compare reasonably well to experimental data. 

With all parameters fixed, integrated and differential cross sections are calculated and 
compared to experimental results [p8|-|32]| . The calculated integrated elastic cross sections 
agree with the experimental values within the numerical and experimental errors for a 
wide range of cm. energies starting at about \fs = 20 GeV up to the Tevatron energy 
y/s = 1800 GeV. Our model does not distinguish between pp and pp scattering or ptt + 
and pn~ scattering, respectively. In the approximation we use here we have C = +1 
exchange only. 

Furthermore the rise of the integrated cross sections in inelastic diffractive reactions as a 
function of \fs is calculated. Our calculated ratio cr s d/(cr e i + cr s d) is in rough agreement 
with experiment. The experimentally observed behaviour that the diffractive dissociation 
part of the cross section increases more slowly with ^Js than the elastic one is reproduced 
qualitatively in our calculation. 

For photo- and electroproduction of p° mesons we restrict ourselves to ratios of integrated 
cross sections. The ratio of single diffractive dissociation to elastic cross sections is found 
to be in the range 0.5 . . . 0.6 for Q 2 < 35 GeV 2 , consistent with the experimental results. 
The Q 2 -dependence of the integrated single diffractive dissociation cross section as well as 
the polarisation ratio of the produced p° mesons also agree with the experimental values. 
To summarise, our model is quite well suited to describe inelastic diffractive hadronic 
reactions at high cm. energies (i/i > 20 GeV) and small momentum transfer. Fur- 
ther progress could be made when including higher cumulant terms in ( |22| ) which would 
contribute to both C = +1 and C — — 1 exchange. 

Another area where our model can be applied is double diffractive dissociation, where 
C = — 1 contributions (odderon exchange) occur already in our present approximation. 
This will be the topic of a following investigation. The upcoming experiments at RHIC 
will be a rich source for new experimental data for both single and double diffractive 
dissociation in hadronic reactions at high cm. energies. Therefore the study of inelastic 
diffractive scattering will remain an interesting and instructive field of work, where effects 
of nonperturbative QCD can be studied. 
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A Definition of 

The light-cone momenta of the quark and the antiquark (or diquark) of the diffractive 
final state X are given by 

P- = Z'P-, P- = (I-Z')P-, (r 7) 

Pt = z'P T + A T , p' T = (l-z')P T -A T , {D,) 

where P_ is the longitudinal and Pt the transverse momentum of X, which moves in 
negative x 3 direction, therefore P + ~ 0. The longitudinal momentum fraction carried by 
the quark is z' . Lorentz invariance requires z' to appear also in the transverse components 
of pt and p' T as defined above. The relative transverse momentum between the quark 
and the antiquark (or diquark) is given by 

A T = + (1 - A P T . (58) 



B Calculation of the form factor 

Starting point for the form factor calculation is the matrix element of the third component 
of the isospin current at x = 

4 = (h 3 (P'M(0)MP)) (59) 

with 

Here Ylip denotes the sum over quark fields u, d and r 3 is the third Pauli isospin matrix. 
The hadrons hi,h 3 are supposed to move in positive x 3 -direction with P + = P' + — >• oo 
(see fK]) . In analogy to the description of hadron-hadron scattering in [|TT|, ^] we therefore 
denote the incoming hadron by h\ and the outgoing hadron by h 3 . The steps required to 
compute the form factor from this expression are completely analogous to those discussed 
in [|11], that lead to the T-matrix element (fj), with the difference being that there are 
now additional contractions between the quarks and diquarks (or antiquarks in the case 
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of mesons) of the hadrons hi, h 3 and the quark fields of the current j% when applying the 
LSZ reduction formalism. By considering the isospin current we ensure that contributions 
which contain subdiagrams arising from contractions between the quark fields of the 
current drop out because they are proportional to Tr r 3 = 0. We describe now the form 
factor calculation for the tt + meson, modeled as ud wave packet. 

by first calculating the matrix 



Using the same notation as in |20| we obtain J 3 ( |59D 
element of j 3 between qq states and then folding with the wave functions of the wave 
packets. 



J 3 



d 2 A 



3T ' ^WW 1 ^ 5 

1 1 



d 2 A 1T 



dzi 



(ud\jg\ud) 



(2 7r )3/2 ^ s bi,-40i{zi, Ait)^5 AiA[ (ud\j£\ud) 
(u(p 3 , s 3 , A 3 )d(p' 3 , 4 A' 3 M(0)\u( Pl , s u A 1 )d(j/ 1 , s[, A[)), (61) 



where Si,Ai are spin and colour indices, respectively and <px j3 are the Fourier transforms 
of the wave functions (|24[) 



(62) 



Applying the LSZ reduction formalism we can express the matrix element (ud\j 3 \ud) (|6l"l) 
as an integral over the quark 6-point-function. We get here only two terms depicted 
graphically in Fig. |ll] which are to be interpreted as follows. We consider a fixed gluon 
background. The quark and antiquark travel in this background and the current either 
hooks onto the quark line (Fig. |TTla) or the antiquark line (Fig. [TT]b) . The matrix ele- 
ment (^T|) is then obtained by averaging over all gluon potentials with a measure given 
by the functional integral (see p0| for the details). In the high energy limit for u und 
d the scattering amplitudes in the fixed gluon background reduce to Wegner- Wilson line 
operators which are closed to a loop W+ by the meson wave functions. This is indicated 



by the dashed lines in Fig. |ll]. Putting everything together we obtain 



Pi+n 



- J^dz J d 2 x T y?3(2:,x r )y3i(2:,x T ) 



el (l- 2 )q T .x T + e -i 2 q T -x T />V + (_ XT)Xt 



(63) 



G 



A straightforward calculation in the MSV gives that the expectation value of the corre- 
lation function of one Wegner- Wilson loop is equal to 1. By a shift in the integration 
variable the d 2 XT integration can then be reduced to a Gaussian integral over the wave 
functions and we find as the final result for the matrix element (|59|) 



(* + (p'Mm* + (p)) 



Pi+< 



dz2z(l 



z) e 



-S 2 a 2 



(64) 
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Figure 11: The two contributions to the matrix element (ud\j^\ud) 



Let us now turn to the proton form factors. In our simple ansatz the proton consists of 
a quark and a scalar diquark, which should be favoured over the vector diquark due to 
dynamical reasons |^9[ . The spin of the proton is then carried by the quark. This together 
with the spin conservation on the parton level has as consequence that, in our model, we 
get for the matrix element of j% between proton states an expression similar to (|64]): 



(p(P' lS ')\M0)\p(P,s)} = P + n%xlF lv {Q 2 



(65) 



with Fi v (Q 2 ) given in ([)!]). Thus we get only a spin non flip and no spin flip contribution 
in the matrix element (f47|), that is, out model gives F2 V (Q 2 ) = 0. This is certainly not a 
very good approximation. But on the other hand the spin flip part in ([|7]) is suppressed 
by \qr\/(2M p ) for — > 0. Thus the matrix element ( |47D is still reasonably described by 
the model for small enough |qr|- 

To summarise, we have outlined in this appendix a calculation of isovector form factors 
using the same methods as for the scattering processes. The results are in essence as 
ro| taking our simple ansatz for the wave functions of the hadrons into account. 



in 
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